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MOND predicts a number of laws that galactic dynamics should obey. I give here a thorough 
description of these laws, their validity domains, and details of their derivation. We do not know 
which of the existing MOND theories, if any, is a step in the right direction towards the inevitable 
deeper MOND theory. It is thus important to pinpoint predictions that follow from only the basic 
premises of MOND: departure from Newtonian dynamics at accelerations a < ao, and space-time 
scale invariance in the limit a <C «o- Such predictions will be shared by all MOND theories that 
■ embody these tenets. The emphasis here is thus on showing how, and to what extent, these MOND 

laws follow from only its basic tenets. It is also important to identify predictions on which MOND 
formulations differ, useful in discriminating between theories. The laws listed here should be obeyed 
£NJ ■ by galactic systems irrespective of their complicated, haphazard, and mostly unknowable histories, 

as Kepler's laws are obeyed by planetary systems irrespective of their complicated formation and 
, evolution. In contradistinction, in the Newtonian-dynamics-plus-dark-matter paradigm, the validity 

^ ■ of such clear-cut laws-which tightly constrain baryons, 'dark matter', and their mutual relations-is 

contrary to expectations. These laws are independent in the sense that no subset of them follow 
from the rest if interpreted within this later paradigm, where they would thus each require a separate 
explanation. Some of these laws hinge on ao in various roles that would seem unrelated if not unified 
by MOND. 
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I. INTRODUCTION 



MOND is an alternative paradigm of dynamics, replacing Newtonian dynamics (ND) plus dark matter (DM) in 
the nonrelativistic (NR) regime, and general relativity (GR) in the relativistic regime [1] (for an extensive, recent 
review of MOND see Ref. [2]). The basic tenets of NR MOND are: (1) A new constant, ao, with the dimensions 
of acceleration, is introduced into dynamics. (2) A MOND theory is required to reduce to standard dynamics in the 
limit of accelerations 1 much larger than ao-formally effected in a given theory by taking the limit ao — )• 0. (3) In the 
opposite, deep-MOND limit (DML), of accelerations everywhere much smaller than ao-formally effected by taking 
the limit ao — > oo, G — > 0, with £1q = Gao fixed-the theory has to become space-time scale invariant (SI), namely 
\Q • invariant under (i, r) — » A(i, r). 

lO ! The original formulation of the DML in [1] posited a relation between the MOND acceleration, g, and the Newtonian 

■ acceleration, g N ~ MG/r 2 , of the form g ~ (ao^iv) 1 , which satisfies tenet (3) (both sides scaling as A" 1 ). However, 
^vq | the third tenet as formulated above is precise and general, and should thus be used instead of the rough relation 
7— ( , between accelerations, which cannot hold exactly, in general. 

We can liken the role of ao in MOND to that of H in quantum mechanics, or to that of c in relativity: All these 
constants, firstly, delineate the validity domain of the old, classical theory, in which they do not appear. Secondly, 

■ they appear saliently in disparate phenomena in the revised domain. 2 MOND tenet (2) is then analogous to the 
correspondence principle in quantum dynamics, for h — > 0, and to relativity tending to Newtonian dynamics for 

r^J c — > oo. 

Beyond the basic tenets, one wishes to construct a NR theory of dynamics based on them, and then extend the 
theory to a replacement of GR. There are several examples of such theories. In the NR regime, we have the suitably 
chosen nonlinear extension of the Poisson equation [3], and a quasilinear MOND formulation (QUMOND) [4]. These 
are classified as 'modified gravity' (MG) theories. There are also 'modified inertia' (MI) formulations, for which there 
isn't yet a full fledged theory, but, on whose implications much can be said nonetheless (see [5, 6]). For more details 
on these and on relativistic MOND theories see the review [2] . 

The departures of quantum theory and of relativity from the classical pictures that preceded them were not mere 
changes in form of the equations of classical dynamics. They each brought totaly new concepts to underlie dynamics. 



* An extended version of a talk given at the 13th Marcel Grossmann meeting, Stockholm 2012 

1 Here, and in many other instances, when I refer to 'accelerations' in comparison with ao, I mean 'all quantities with the dimensions of 
acceleration'. This is analogous to relativity reducing to Newtonian dynamics when all quantities with the dimensions of velocity are 
much smaller than the speed of light (not only velocities, but also, e.g, the square root of the gravitational potential). 

2 As we shall see below, it is not ao itself that appears in deep-MOND phenomena: it is always the constant Oq = Gao that does. 
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This may well turn out to be the case for an eventual deeper MOND theory. A hint of such eventuality may be seen 
in the possible connections that MOND brings out between local dynamics and the global state of the universe (see 
Ref. [7] and references therein). 

Existing MOND theories are extensions of the pre-MOND dynamics (NR and relativistic) that introduce ao, add 
new degrees of freedom, and modify the underlying action; but, they do not deeply depart in spirit from their 
predecessors. 3 Perhaps one of them will turn out to be an effective theory that captures the essence of the deeper 
MOND theory. 4 But perhaps none of the existing MOND theories will turn out to be a step in the right direction; 
perhaps their merits are summed up in their being self-consistent theories that embody the basic tenets of MOND. 

It is thus noteworthy that one can predict a number of general laws, based almost entirely on the basic premises. 
These laws are important to recognize and understand in detail for several reasons: (a) They constitute the core 
predictions of MOND, obeyed by all MOND theories, present and future, that embody the basic premises, (b) They 
focus attention on well-defined, easy-to-grasp, and memorable sub-predictions of MOND, subsumed by the general 
prediction that MOND would make of the general dynamics of a system, (c) They involve only limited information 
about the systems (such as global attributes) and hence are easier to test on large samples [in contrast with, e.g., 
rotation-curve (RC) analysis, which requires knowledge of the full RC, and the detailed baryon mass distribution], 
(d) They constitute a 'to-derive' list for the competing paradigm of ND-plus-DM. (e) Pinpointing and quantifying 
the remaining differences in these laws as predicted by different theories will help discern between theories. 

Succinctly formulated, the laws recognized to date arc: 

(i) Speeds along an orbit around any bounded mass, M, become asymptotically independent of the size of the orbit 
(being dependent only the shape of the orbit). For example, for circular orbits V(r — > oo) —¥ V 00 (M). This 
replaces Kepler's 3rd law T 2 = K(M)r 3 with T = K(M)r. 

(ii) V X {M) = (MGao) 1 / 4 , which dictates the dependence of K on the central mass. 

(iii) A mass discrepancy appears when we cross a = a . In disc galaxies, the transition from 'baryon dominance' to 
'DM dominance' occurs always around the radius where V 2 (r)/r = a - 

(iv) Quasi-isothermal systems, which well model many galactic systems, have mean surface densities XI ^ S M = 
a /2irG. 

(v) Quasi-isothermal, or deep-MOND, systems of mass M, have characteristic velocity dispersion a ~ (MGao) 1 ^ 4 . 

(vi) The external field in which a system is falling affects its intrinsic dynamics: the external-field effect. 

(vii) Disc galaxies behave as if they have both disc and spherical 'DM' components of predictable properties. 

(viii) MOND endows self gravitating systems with an increased, but limited stability. 

(ix) The incremental acceleration MOND predicts (that attributed to 'DM') can never much exceed ao- 

(x) The central surface density of 'dark halos' is < S M . 

It is reassuring that even the more elaborate prediction of RCs in MOND follows almost entirely from the basic 
premisses alone. 

The above-listed laws have been discussed before [2, 7, 8]. But, in the following sections I consider in great detail 
their meaning, implications, and how they are derived. And, in particular, I will show to what extent, and how, 
these laws follow from the basic tenets alone, and what aspects of them are theory dependent or require additional 
assumptions. 

Some generalities to take home concerning these laws are: 

a. To a large extent they do follow from only the basic tenets, and are thus shared in one way or another, by all 
MOND theories that embody these tenets. 

b. Several of them revolve around a in different roles. 

c. These laws are independent as phenomenological laws. 5 For example, if interpreted as effects of DM in ND, one 



3 Such theories might be analogous to the pre GR attempts to relativitize gravity by writing Lorentz invariant theories for the gravitational 
potential— as, e..g., in Nordstrom's theories— before the full force of gravity-as-space-time-geometry has been appreciated. 

4 The appearance in all these theories of an interpolating function is an indication that, at best, they can be effective theories. 

5 In the same sense that without the unifying framework of quantum dynamics, the different quantum phenomena-such as the black- 
body spectrum, the photoelectric effect, the hydrogen-atom spectrum, superconductivity, etc.-would appear unrelated phenomena that 
somehow involve the same constant h. 
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can construct model families of baryons plus DM that will satisfy any subset of these laws but not the others. 

d. Some of these laws, when interpreted in terms of DM, would describe properties of the "DM" alone [e.g., (i), (ix), 

and (x)], of the baryons alone [e.g., (iv)], or relations between the two [e.g., (ii), (iii), and (vii)]. 

Comparisons of these predicted laws with data are discussed at length in [2] and in many previous studies, and will 
not be considered here. 



II. ADDITIONAL ASSUMPTIONS 

In deriving the MOND laws, I assume that MOND does not involve additional dimensioned constants. 6 I also 
assume that it does not involve dimensionless constants that much differ from unity, fc, whose appearance is 
tantamount to having in the theory several, greatly different, acceleration constants, ao, kao, k ao, etc., that may 
play different phenomenological roles. My whole discussion below hinges on the premise that this does not happen. 7 
This assumption means, in particular, that the DML (with its scale invariance) applies already below a relatively 
large fraction of ao- It also implies that Newtonian dynamics holds approximately already above a few ao- I will 
call this latter 'weak compatibility' of MOND with standard dynamics. 8 In some special cases I shall assume 'strong 
compatibility': the requirement that MOND approaches standard dynamics very fast, beyond ao; more quantitatively, 
that for g 3> a , |g — g N \/g = o(a /g), namely, that |g — g N \/a — ► for g/a — > oo. 9 Such strong compatibility 
is hardly ever relevant to galactic dynamics, which only probe accelerations up to g ~ lOao- However, dynamics in 
systems such as the solar-system and binary pulsars, do probe accurately much higher accelerations. 

Existing MOND theories all involve some interpolating function between the Newtonian and MOND regimes, of 
the form fi(A/ao), where A is a quantity with the dimensions of acceleration. Weak compatibility means, for these 
theories, that (i(x) ~ 1 already beyond x ~ a few. Strong compatibility means that [1 — fJ,(x)} <C p/x for x ^> I (p is 
some constant of order unity). 

More assumptions are usually made implicitly, and I spell them out here. It is assumed, for example, that, as 
in Newtonian dynamics, all the (NR) effects of a body at distances much larger than its size depend only on the 
mass of the body, not, e.g., on its internal structure, or on attributes other than its total mass. 10 This means, for 
example, that predicting the dynamics of a system made of bodies small compared with their separations requires 
only knowledge of their masses. I shall refer to such bodies in such a context as pointlike. This assumption can be 
shown to hold in the known NR MOND theories, such as the nonlinear Poisson theory [3], QUMOND [4], and the 
NR limits of all the relativistic MOND theories discussed to date. But, it might be that a future MOND theory 
will violate this assumption (e.g., by adding coupling constants other than the mass). This requires , in particular, 
that the bodies and systems we are dealing with have non-zero masses. There are theories and situations, where the 
effective gravitating mass of a system might vanish, as happens, e.g., in the NR limit of BIMOND, for a system of 
equal amounts of matter and twin matter (see Ref. [13] for a detailed discussion). In this case, which requires special 
treatment, and which I bar in this paper, the asymptotic field does depend on the structure of the system even if it 
is very small in extent (see below). 

III. MORE ON THE DEEP-MOND LIMIT 

The DML is pivotal in understanding the MOND laws; so I expand a little on its properties. The dynamics in this 
limit was originally [1] epitomized by the relation g ~ (aogw) 1 ^ 2 between the MOND and Newtonian accelerations of 



6 Counter examples are the MOND theories discussed in Ref. [9], which adds a length scale, and in Ref. [10], based on a suggestion in 
Ref. [11], which adds a velocity constant. 

7 TeVeS [12], and its NR limit, is a counterexample: it does have additional dimensionless parameter(s). TeVeS does not tend to GR 
in the mere limit ao — > 0, and to make it agree with various high-acceleration constrains, one of these constants is forced to be much 
different from unity. 

8 The second tenet constitutes a compatibility requirement, but it requires compatibility only in the very limit ao — > 0. 

9 This definition of strong compatibility departs from the strict generality I am adopting everywhere else. We are not supposed to use 
accelerations in defining our assumptions, since it is not clear that a MOND theory, and the criteria for its Newtonian compatibility, 
involves (only) accelerations. Even in existing theories, g is not a function of g N , which this formulation of the criterion assumes. For 
example, in MI theories, we may have Newtonian behavior for some orbits, and non-Newtonian behavior for others, even if they all 
have similar accelerations. So in the rare cases where I will assume strong compatibility it should be understood to apply only to the 
problem at hand (e.g., to the types of orbits under discussion). 

Namely, finite size effects are assumed to decay with distance faster than the dominant contribution that depends only on the total 
mass. 
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test particles. However, except for special circumstances, 11 the MOND acceleration cannot simply be a function of 
the local Newtonian acceleration. Such an algebraic relation does not hold exactly in any known MOND theory, and 
is also not consistent with conservation laws. In MI theories the MOND acceleration, unlike g w , is not even a function 
of position, but depends also on the full trajectory. So, unlike the formulation of the DML here, by the requirement 
of scale invariance, which is exact and clear-cut in meaning, the g — g N relation is neither. Still, this algebraic relation 
is very useful, and does apply under some circumstances, where it can be derived from the above basic assumptions, 
as we shall see. 

A theory that is space-time scale invariant can be written so that all its dimensioned constants do not change 
under scaling of the length and time units by the same factor. Specifically, in our case, the DML is invariant under 
t — > Xt, I — » XI (applied only to the theory's degrees of freedom, not to the constants G and a ). In addition, like any 
theory, the DML is invariant to the change of units (which changes also the constants): t — >• at, I — > ai, under which 
ao — > a ao, G — > aG. Taking A = a -1 , and applying both, tells us that the theory is invariant under ao — > a ao, 
G — > aG. This means that only the product 17 o = Gao can appear in the theory: 12 In a universe governed strictly 
by a DML theory, there is neither ao nor G, only 17o appears. If it were not for the existence of a Newtonian regime 
of observed phenomena, we would not have known of G, and would not have spoken of ao. Below, I shall use 17o in 
deep-MOND results. 13 It is also enlightening to write DML results in terms of the MOND mass M = c 4 /S7 .[7, 8] 

Now that we know that only 17o appears in the dynamics, we infer (e.g., [7]) that the theory is invariant, more 
generally, to all scalings that do not change the value of Slo (when considered as change of units): This is a two- 
parameter family I — > a£, t — > $~ x t, m — > (af3) 4 m. This specific two-parameter invariance comes about because NR 
MOND involves only the constant ao in addition to G, together with its scale invariance. 14 Thus, if rj(i) is a system 
history for m i; then aVi{f3t), with velocities a/3vi(j3t), is a system history for masses {af3) A m,i (with appropriately 
scaled initial conditions). 15 

Imagine then an isolated DML system, in virial equilibrium, of total mass M made of test particles of masses 
m, = Mqi, and of characteristic size R, and characteristic internal speed V; so we can write for the particle speeds 
Vj = VuJi, with some average of the cot (say their mass weighted rms) being 1. Define 

Q = V2 = V2 I R = I (!) 
(M17Q) 1 / 2 (aoMG/i? 2 ) 1 / 2 (a Q g N y/ 2 ' 1 ' 

where g = V 2 /R is what we would identify with the characteristic MOND acceleration, and g N = MG/R 2 is the 
characteristic Newtonian acceleration. It follows from the above that every such system is a member of a two-parameter 
family of systems having all values of M and R (provided they are still in the DML; i.e., with MG/R 2 <C a ), but 
all having the same value of Q. In other words, the value of Q for a DML system cannot depend on any of the 
dimensioned attributes of the system (M, R, V), only on dimensionless and scale-invariant attributes, 16 such as g,, 
oji, shape parameters, orbit distribution, etc. all of which are the same within each family of systems. For example, 
DML, pure-exponential discs of all masses and sizes are in the same family; all DML isothermal spheres (ISs) with 
the same constant anisotropy ratio are in the same family, etc. 

While Q is predicted to be strictly constant within each family, the above arguments alone do not prevent it from 
varying greatly among the families. Indeed, many of the dimensionless parameters on which Q is permitted to depend 
can take up values much smaller or much larger than unity, and very different among the families. Still, I argue that 



For example, for one dimensional configurations in some theories. 

The theory's degrees of freedom can be normalized so that they transform under scaling according to their dimensions, and I adopt this 
normalization. For example, the DML of the nonlinear Poisson formulation of MOND can be written as a = — V0, with V ■ [| V0|V</>] = 
4nGaop, where cj> has dimensions of £ 2 /t 2 , and indeed remains intact under scaling; this is the desired normalization. But we could 
also work with ip = (p/ao, with the equations taking the equivalent form a = —aoS/tfi, with V • [|V?/>|Vi/>] = (4irG/ao)p. Now, ip has 
dimensions of i but must still remain intact under scaling to retain scale invariance, not consistent with its dimensions. Indeed we see 
that with this normalization, ao and G do appear separately, not as Gao. In MI theories, we have schematically, aa/ao = — V</>, with 
A0 = AnGp. To have scale invariance, <fi has to scale as A -1 , which is incongruent with its dimensions-and indeed G and ao appear 
separately. Working with tp = ao<f>, for which the scaling dimension is the desired one, we have aa = — V?/> with Alp = AwGaop. 
In Newtonian dynamics, G 1//2 has to be introduced as a dimensioned proportionality factor between the inertial (m) and the gravitational 
mass (mc). The former is defined so that inertial forces are ma, and the latter is defined such that gravitational forces scale as m^/r 2 ; so 
tug has the dimensions of G x ^ 2 m. (The masses, in, we use all along are the inertial ones.) In the DML, gravitational forces must scale as 
p(m, ao, G)/r to retain scale invariance. Dimensionally, we must have p ~ flj/ 2 m s ^ 2 ; so gravitational forces arc then F ~nl /2 m 3 / 2 /r. 
Had we built MOND, for example, around a new length constant, I, instead of an acceleration, and require scale invariance in the limit 
of large distances, I — > 0, then only the ratio G/i would have appeared in the limit. In this case the two-parameter family of invarianccs 
would have had m — ¥ (a/3) 2 m. 

For continuum systems, if the density and velocity fields p(r,t), v(r,t) arc a solution, so arc a/3 4 p(r/a, 0t), a/3v(r/a, /9t). 
So, Q cannot depend on the dimensionless, but non-scale-invariant MG/R 2 ao. 
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due to our extra assumptions, Q cannot vary greatly among sub-ao systems, and is, in fact, of order unity for all of 
these. 

The argument hinges strongly on my above-stated assumption that MOND does not involve dimensionless constants 
very different from unity. 17 In particular, it follows that the transition from the DML to the almost Newtonian 
regime is relatively narrow, and always occurs within a factor of a few do in acceleration. Take then a system with 
MG/R 2 = do- The above assumption says that its dynamics is not much different from Newtonian, so, with our 
definition of V, it satisfies approximately the Newtonian virial relation MG/R <~ V 2 . Eliminating R gives Q ~ 1 for 
such systems. But, for each DML family of scaled systems-which fills the region (much) below the borderline parabola 
M = (a /G)R 2 in the R — M plane-there are members near the borderline. By the above assumption, the dynamics 
of these members is not very different from that of their neighbors on the borderline that have a similar geometry, 
mass distribution, and other dimensionless attributes. 18 So they too satisfy Q ~ 1, and so do all the DML systems. 19 
We have thus derived the approximate DML relation g ~ (aog,^) 1 / 2 between the characteristic global parameters of a 
system from the basic tenets and our additional assumptions. Such a relation does not necessarily apply locally. 

The meaning and the expected dynamics of the relativistic DML are not understood yet: Because it turns out that 
ao w cHo/2tt (Ho is the Hubble constant) [1], a system that is both relativistic, and of accelerations a <C ao has to 
have a size larger than the Hubble distance; so this double limit does not apply to any know system. 



A. Nonlinearity of the deep-MOND limit 

All existing MOND theories have a nonlinear nonrelativistic DML. 20 Is this forced by the basic tenets? The global 
relation g ~ (oo^n) 1 ^ 2 , which we derived from our assumptions, and, e.g., law (ii), which follows from the basic 
tenets, imply already that we cannot superpose accelerations produced by several mass distributions, as is possible in 
Newtonian dynamics (DML accelerations scale as the square root of the mass). But this does not necessarily imply 
that the theory is nonlinear. A linear NR dynamics is one where the motion of a test particle induced by a mass 
distribution p\ + pi can be gotten by superposing the motions calculated for p\ and p2 separately. This means, more 
precisely, that the particle equations of motion of all particles are of the form 

(OrO(i) = (Cp)( ri ), (2) 

where O is a linear (differential, possibly time-nonlocal) operator acting on the trajectory to give another vector 
function of time, (0rj)(t), and C is a vector, linear operator acting on the field p. Only the constants G and ao may 
appear in the two operators. However, since p appears linearly and only on the right hand side, G can only appear 
as its prefactor. Also, for a DML theory we can only have fio appearing. This means that C = OqC, where C, like 
O, cannot contain any dimensioned constants. In particular, since we do not have constants with dimension of time, 
O must be cx d n /dt n . (Time cannot appear explicitly because of time-translation invariance.) Assume, for simplicity, 
that C does not involve time derivatives. Since fioP has dimensions £t~ 4 , C can have no length dimensions, and n = 4. 
(If we do allow time dimension — p for C, with a pth time derivative appearing therein, we must have n = 4 +p.) We 
then must have from the basic tenets 

^ - ^oCp(r), (3) 

where C is dimensionless. This is schematically tantamount to replacing the standard, nonlinear DML relation 
a 2 ~ iloM/r 2 by the linear a ~ £loM/r 3 , which retains scale invariance. 
As one example of many, we could take 



Demensionless system parameters can, of course, take up small or large values. 

For example, all DML, pure-exponential discs have a similar Q value to borderline exponential discs; all DML ISs with the same constant 
anisotropy ratio have their borderline counterparts, etc. 

It is easy to see that if, for example, contrary to our assumption, the DML (i.e., the region of scale invariance) applies only up to 
accelerations kao (k -C 1), above which the theory transits quickly to Newtonian behavior, then DML systems have Q ~ fe; and, in 
principle, fc could even depend on system type. 

Relativistic versions are nonlinear also in the high-acceleration regime, as GR is. 



6 



with 

1 f p(v')d 3 v' 

(The field ip is not a modified gravitational potential, since is not the acceleration.) 

Note that in a theory of this type both the kinetic and gravity terms are modified so it is neither a MG theory, nor 
a pure MI theory. 21 

Equation (3) is the most general equation of motion describing a linear DML theory (for p = 0). It is, however, not 
an acceptable DML theory: 1. It is a high order theory, with all the known associated problems: the awkwardness 
of having to dictate 4 initial conditions, the inevitable presence of instabilities as implied by Ostrogradski's theorem 
[14], discussed in detail in Ref. [15], and the related problem of exploding solutions. 2. There is no general MOND 
theory that has this equation of motion as its DML, and reduces to Newtonian dynamics in the limit oo — > 0. Clearly 
there is no kinetic action that interpolates between {d 2 Ti/ dt 2 ) 2 and the Newtonian {dxi/dt) 2 and involve only ao- It 
was shown generally in Ref. [5] that there is no local kinetic Lagrangian (one depending only on time derivatives of 
Yi) 1 with Galilei invariance, that has both a correct DML and a Newtonian limit. 

Thus, if we exclude theories of this type, nonlinearity is forced by the basic tenets. Still, such a DML formulation 
may be useful as a toy testbed in various contexts. 



IV. MOND LAWS OF GALAXY DYNAMICS 



I now discuss the MOND laws in detail: how they are derived, and to what degree they follow from the basic tenets. 

Some of these laws can be described in terms of the properties of a putative DM distribution that reproduces 
MOND effects. However, not all MOND theories have an equivalent description in terms of DM. In theories that 
only modify the gravitational potential from the Newtonian tp N to the MOND potential <f>, the modification is fully 
described by the (Newtonian) gravitational effects of 'phantom matter' (PM) 22 of density 

Pp = i ^(A0-A^). (6) 

But, this is not the case, for example, in MI theories where the modification may depend not only on position, but 
also on the particle orbits. 

Still, when analyzing the dynamics of galactic systems the data has always been limited, never relying on studies of 
all orbit types, at all positions. With such limited information on the dynamics of the system, a description in terms 
of DM equivalent may be possible even when it does not apply to the full dynamics. I refer here to the properties 
of 'DM' obtained in this limited way. For example, rotation-curve analysis pertains only to the dynamics of circular 
orbits in the equatorial plane of an axisymmetric disc galaxy. With this limited data, the departure from Newtonian 
dynamics can be attributed to a spherical halo of DM, even in MOND theories where this would not be possible for 
general orbits. 



A. Asymptotic constancy of the orbital velocity around an isolated mass 

In MOND, the orbital speeds become independent of the size of the orbit for an orbit far outside some isolated 
mass. This follows simply from scale invariance of the DML: Under scaling the size of the orbit and all the relevant 
times are multiplied by the same factor; so velocities do not change (the scaling of the central mass distribution is 
immaterial for asymptotically far orbits). 

More formally, consider some isolated mass distribution, p(r), of total mass M, having a characteristic radius R 
(some minimal radius containing most of the mass in the system). Take r(t) to be some asymptotic, test-particle 
orbit. By 'asymptotic' in the MOND context one understands two conditions: 1. |r(i)| ^> R; so the orbit is far outside 
the mass, which can then be considered a point mass; so, the dynamics of such orbits are oblivious to the exact form 



The equation of motion (3), with the choice (4), can be gotten from a Lagrangian C <x ~}2i[(d 2 ri/dt 2 ) 2 — 2QoV'( r i)]- 
This is a very useful concept in the context of MOND. For some applications see Refs. [16-19]. 
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of p or its extent R. 2. |r(t)| 3> r M , where r M = (MG/ao) 1 ^ 2 is the MOND radius of the system; 23 so the orbit is 
wholly in the DML. In particular, the orbital dynamics does not change if we expand the system so that R^> r M (but 
still |r(i)| ^> R). This puts the whole system is in the DML, and we can apply scale invariance: For any asymptotic 
orbit r(i), the orbit r(i) = Ar(A~ 1 t) is an orbit for the distribution p = A~ 3 ,o(A~ 3 r), also having total mass M. But, 
since the structure of p is immaterial, we conclude that f is an orbit for p(r) itself. The velocities on the scaled orbit 
are v(t) = v(X~ 1 t). Thus, for such asymptotic orbits the orbital velocities are independent of the size of the orbit. 

This is an obvious result for logarithmic potentials, but I have shown that it follows, more generally from the basic 
tenets. For circular orbits, this means that for r 3> i?, r M , the rotational speed becomes independent of r, 24 a result 
that is exact for isolated systems in all MOND theories. 25 

For comparison of this prediction with data see the many published compilations of RCs, summarized in Fig. 15 of 
Rcf. [2] (see also their Figs. 21-27). 

If interpreted in terms of DM, this prediction pertains to a pure-halo property, since asymptotics is governed by 
the strongly dominant putative DM halo. 



B. The mass-asymptotic-velocity relation 

Given that the asymptotic dynamics of test particles depend only on the total mass, and given that the only quantity 
with dimensions of velocity that can be constructed from mass, length, and Vt is (Mflo) 1 ^ 4 , the orbital speed on any 
asymptotic trajectory has to be of the form 

v(t) = (Mn )^t(t) = c C(t), (7) 

where ((t) covers all dimensionless orbits, which are independent of the specific MOND theory, of details of the mass 
distribution, or of the size of the orbit (they are the dimensionless classical orbits in a logarithmic potential). In 
particular, for circular (constant speed) orbits, \C(t) \ is a universal constant for the theory. The common convention is 
to normalize fio (or ao, given G), such that for circular orbits \C(t)\ = 1. Thus, the mass-asymptotic-velocity relation 

(V \ 1/4 

M = O "Vi = M i^-fj , (8) 

between the total central mass and the asymptotic circular velocity, is an exact prediction of the basic tenets of 
MOND, and of all the theories that embody them. 

For recent comparisons with data see, e.g., Refs. [20, 21], and Figs. 3-4 in Ref. [2]. 

If interpreted in the context of DM, this prediction constitutes a relation between the total baryonic mass and a 
halo property, the quasi-asymptotic speed. 



C. The mass discrepancy correlates with acceleration, and in disc galaxies appears always at a fixed 

acceleration value 

In disc galaxies for which a = V 2 (r)/r > ao in the inner parts, the mass discrepancy is predicted by the basic tenets 
to appears always around the radius where a — ao- Because the only quantity with the dimensions of acceleration 



In systems with R < r M -namety, when most of the mass is contained well within Tm~ the MOND radius plays the role of a transition 
radius between a Newtonian space region r < r M and the MOND space region r > r M - Interesting change of appearance may then be 
observed across this radius. When R 2> r M there is no such transition, as the whole system is then in the MOND regime. This aspect 
of Tm is similar to that of the Schwarzschild radius, r s = 2MG/c 2 , in relativity, where for R 3> r s the system is wholly in the classical 
regime, but when R < r s , r s does mark a transition radius. 

I assume, that the theory is not so odd as to permit orbits with the same radius, but different velocities. 

Interesting departures occur when the theory allows gravitating masses of opposite signs, as e.g., in BIMOND with twin matter [13]. In 
a system with vanishing total mass, the asymptotic speeds are not constant: Here, the scaling of p(r) itself does affect the asymptotic 
field (the asymptotic field is no more the canonical MOND logarithmic potential). For example, in NR BIMOND, the field equation was 
solved exactly in Ref. [13] for two opposite masses ±m in the DML . Asymptotically, the potential is <j> -(mfioj'^r ■ d/r 2 , where 
d is the dipole separation. So, asymptotic speeds decrease as (mC2o) 1 / 4 (c(/r) 1 / 2 . The breakdown of the general result occurs because 
the asymptotic potential is not invariant to the scaling of the mass distribution (under which to d -> Ad). However, masses of opposite 
signs repel each other, so one cannot have a self gravitating system of this type. In what follows I shall ignore such systems. 
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in this context is V 2 (r)/r, the transition from Newtonian to MONDian dynamics, dictated by the basic tenets, must 
occur at a radius where V 2 (r)/r = ao, if such a radius exist. Systems in which V 2 /r < cio everywhere are predicted 
to show MONDian behavior everywhere, with no transition. A similar behavior is expected in existing theories also 
for quasi-spherical, pressure-supported systems. 

For comparison with data see, e.g., Fig. 12 in Ref. [22] for pressure-supported systems, and Fig. 10 of Ref. [2] for 
disc galaxies, both showing clearly that the mass discrepancy is a function of a/ao with (gradual) onset at a = ao. 

If interpreted within the DM paradigm, this prediction constitutes a tight connection between visible matter and 
DM, since it predicts where DM dominance takes over from visible-matter dominance in a any system. 



D. Quasi-isothermal spheres have mean surface densities (£) < S M = ao/2-wG 

Self-gravitating ISs of finite mass (possibly with anisotropic velocity distributions) do not exist in Newtonian 
dynamics: If the mass of a self gravitating body is finite, then beyond a certain radius its enclosed mass is nearly 
constant. Beyond this radius, Newtonian speeds of bound particles have to decrease as r -1 / 2 , which is inconsistent 
with isothermality. 26 In MOND, however, asymptotic isothermality for a finite mass is natural, as we saw in section 
IV A. It follows that for a self gravitating sphere to exist and be protected from evaporation (by MOND), it cannot 
be deep in the high-acceleration regime; namely, it cannot have a size R <C r M . But R > r M means that the mean 
surface density (E) - M/ttR 2 < E M . 27 

The above results follow from only the basic tenets of MOND, but Em should be taken as only a rough upper 
limit on, and an accumulation value of, (E). MOND ISs have been considered in detail for the nonlinear Poisson 
formulation in [23]. (These results apply also in QUMOND, since the two theories coincide for spherical systems.) It 
was shown there (Fig. 2) that IS families with different anisotropy ratios have upper limits (E) < aE M , with a ~ 1 for 
isotropic orbits, becoming somewhat larger for more tangential orbit, and becoming smaller, even as small asa = 0.1 
or smaller for more radial orbits. 

Inasmuch as globular clusters, dwarf spheroidal galaxies, elliptical galaxies, and galaxy clusters are quasi-isothermal 
(not all individual systems within each class are approximately IS), we can apply this prediction to them. Indeed they 
satisfy this inequality, as is shown, e.g., by the 'Fish law' for ellipticals, and see Fig. 6 of Ref. [22]. This prediction 
might even be pertinent to giant molecular clouds [24]. 

This law, if interpreted within the DM paradigm pertains to pure baryonic attributes. 



E. A mass-bulk-velocity relation 

For isolated, virialized, self-gravitating systems with (S) < E M , the basic tenets of MOND predict a correlation: 

" i ~ Mn ° =<='(£>■ (9) 

where M is the total mass, and a is the characteristic, bulk velocity. 28 This follows from the results of section III 
identifying a with V there. 

Despite the similar appearance, this law is very different from law (ii), — MQ,q: (a) The latter is an exact 
relation, the former only a correlation with possibly large scatter, (b) The latter involves the asymptotic, circular 
speed, the former a bulk, characteristic speed within the system, (c) The latter applies to all isolated bodies, the 
former is predicted for limited classes. 29 

Quasi-isothermal spheres are a relevant special case, since by law (iv) they have (E) < E M -they are either DML or 
borderline cases and so arc all predicted to satisfy correlation (9). Seen differently, in ISs the bulk a is, by definition, 
also the asymptotic a. But the latter must be ~ for the system. Thus, in this case the M — o correlation does 
follow from law (ii). 



Viewed differently, for an IS, the Newtonian equation of equilibrium is p = dln(p) / dln(r) oc —M(r)/r. Beyond the radius containing 
most of the mass, p thus tends to zero, while finitcness of the mass would imply that asymptotically p < —3. 
Non-isothermal spheres, such as stars, can, of course, have arbitrarily high values of (E). 
Defined, for example, as the mass-weighted RMS velocity a 2 = (V 2 ) = M~ 1 J p(r)v 2 (r)). 

For example, for an isolated Sun-like star, which is highly Newtonian in its bulk ((E) ^> Em), law (ii) should still hold, with Voo ~ 
0.4km s _1 . However, the characteristic bulk velocity is a f» (r M / R) 1 ^ Voa ~ 400km s _1 , so clearly the mass-bulk-velocity relation is 
not satisfied. 



9 



Clearly, there is considerable range that is allowed for the value of Q — a 2 /(MQo) 1 / 2 by the arguments leading 
to correlation (9). But this range in the ratio is small compared with the relative spread of values of a 2 or M 1 / 2 
separately within the gamut of systems to which this prediction applies-from globular clusters and dwarf spheroidals, 
through galaxies of all types, to clusters of galaxies. 

In specific MOND formulations we may be able to derive more concrete versions of this prediction. Some such 
predictions are known. In the DML of modified- Poisson theory of Bekenstein and Milgrom [3], and of QUMOND [4], 
all isolated virial systems satisfy 

2 = ^1-E^ 2 )' (10) 

where m = rrii/M are ratios of constituent masses, m, to M, and a is the mass- weighted root-mean-squared (3-D) 
speed. This, more concrete, and exact relation follows from some additional symmetry of the DML of these theories, 
which is not part of the basic tenets [4, 25, 26]. For a smooth mass distribution (made of 'test particles', for which 

Qi^ 2 ^ 1) wc g°t a universal value Q = 2/3 for any such DML system irrespective of its structure. 

We do not have such a general result for MI theories. However, for such theories there is a simple prescription for 
calculating the rotation curve of any axisymmetric disc galaxy [5]. Taking advantage of this, I calculated [27] the Q 
value for large classes of thin DML disc galaxies, where only the rotational speeds contribute to a (defined as above). 
I found, remarkably, that all have Q values within a narrow range: Q « 0.73 ± 0.01. 

Isothermal spheres with a constant but arbitrary velocity-anisotropy ratio were considered in Ref. [23] for the full 
MOND (not only DML), nonlinear Poisson theory (the results apply also to the later QUMOND). It was found that 
for all such spheres 2/3 < Q < 1, where the exact value depends on the anisotropy ratio, and on the closeness to the 
DML. 

In the toy, linear DML theory described by eq. (3) with the choice (4)-which satisfies the third tenet-one can easily 
derive an exact relation of the form 

(a 2 HtW-i(£-^) t , (11) 

i<j ' 1 

for a steady-state system, where (a 2 ) is the time average, mass-weighted system average of the squared acceleration, 
and ()t signifies time average. If we replace the sum with ~ M 2 /R 2 , and (a 2 ) with ~ <j 4 /R 2 , we get Q <~ 1 for this 
'theory'. 

This MOND law underlies the observed Faber- Jackson relation for elliptical galaxies, and its generalization to all 
pressure-supported systems, as discussed by Sanders [28]. He applied it to systems spanning almost 10 orders in M and 
a 4 , and also extended the discussion to a three-parameter fundamental plane, allowing for an additional parameter 
to M and a. For additional comparison with data see Ref. [2]. In addition to the intrinsic scatter expected around 
relation (9), there is scatter due to conversion of line-of-sight dispersions to 3-D ones, to other deviations from the 
definition of a used in our derivation, to departures from assumptions in the derivation (e.g., isothermality) , to errors 
in the masses (which have to be converted from luminosities), etc. 

When interpreted in the DM paradigm this prediction pertains sometimes to pure-baryon properties (when (£) ~ 
£ M , so baryons dominate), and sometimes to a relation between baryons and 'DM' (when (£) <C S M , so 'DM' 
dominates). 

F. The external-field effect 

Unlike Newtonian dynamics, MOND is nonlinear even in the nonrelativistic regime. It generally does not satisfy the 
strong equivalent principle; so effects of an overall acceleration on the internal dynamics of a system are gencrically 
expected. To be able to say what constraints the basic tenets impose on such effects I have to confine myself here to 
theories whereby only the instantaneous value of the external acceleration matters. This excludes from the discussion 
a large class of MI theories that are time nonlocal [5]. In these, the full (external) trajectory of the system enters, 
which complicates the discussion. Some of the possible consequences of such nonlocality are discussed briefly in Ref. 
[6], but what follows here does not apply to such theories. 

Consider then a system of mass m ('system to'), and extent r, that is falling in the field of a mother system with 
acceleration whose instantaneous value is g . Assume that the theory and conditions are such that, to a good enough 
approximation, all the information about the mother system enters the dynamics within m only through g . One can 
then write 



a = a(m, r, a , G, g , n , a), 



(12) 
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where a stands for the internal acceleration runs of elements of m, namely the full acceleration in the field of the 
mother system minus g (suppressing the dependence on position, time, and particle index). It is written as a function 
of all the available dimensioned independent parameters, as well as of n , the unit vector in the direction of g , and 
of a, which stands for the many dimensionless parameters that characterize the configuration, such as all the mass 
ratios, and all the geometrical parameters (angles, ratios of all distances to r, etc.). Here I am only interested in 
scaling laws of the dimensioned parameters-for example, in how |a| depends on the dimensioned system attributes -so 
I shall suppress the dependence on n and a. 

Since a/g is dimensionless, it can depend only on dimensionless quantities; so we can write, most generally 

a = g F(i],9), > 9=—. (13) 

r 2 9o 9o ao 

When go <C |a|, its effects can be neglected. So here I shall be interested in the opposite case, of external-acceleration 
dominance, g » |a|. 30 The above choice of dimensionless variables is useful for this case. Clearly, F*(0, 9) = 0. So 
we are interested in the behavior of F* to lowest order in 77. We shall see below that external-field dominance requires 
77 <C 1 when 9 > 1, and the scale- invariant condition n <C when 9 < 1 (in which case the whole problem is in the 
DML; i] and 9 both scale like A -1 under scaling); so we can write this condition generally as 77 <C min(l, 9). 

We do not know that a MOND theory is necessarily expandable in powers of 77 near 77 = 0. But assuming that it 
does, I write 

a « g^H (9), 77<mm(l,0). (14) 

(I assume that q does not depend on 9; see below.) 

To constrain q and f (9) I now employ the basic tenets of MOND. The limit a — > 0, namely when ao <C |a| <C go, 
is strongly Newtonian for all accelerations, and is within the validity domain of eq.(14). For this region, c/o 

and ao 

have to disappear from expression (14). This implies that q = 1, and that |f(0 3> 1)| ~ 1, such that (mG/r 2 )f (00) 
is the internal Newtonian acceleration field within m. This means that the internal dynamics is Newtonian for any 
value of g N when 0^1; i.e., also when g N <C a . In other words: whenever the external field is highly Newtonian 
and dominates over the internal field, the latter is necessarily Newtonian. This result holds also when q depends on 
6, because then we still must have q(9 — > 00) — > 1. 

More generally, inasmuch as q = 1 for all 9, we can write eq.(14) in its full validity domain (external-field dominance) 

as 

a=^f(0). (15) 

This means that when the external field is dominant, the internal dynamics is always quasi Newtonian, in the sense 
that the accelerations scale as mG/r 2 , only with an enhanced effective constant G e // ~ G|f(0)|, and with not-quite- 
Newtonian geometrical aspects that stem from the fact that f has different geometric properties than f(oo): for 
example, f depends on the direction relative to n , and on the theory at hand, while f(oo) does not. 

When 9 <C 1 the whole system is in the DMR, where the basic tenets dictate that eq.(15) becomes scale invariant. 
Under scaling, 9 scales like go, namely 9 — > A _1 (since g is a DML acceleration of the mother system it scales as 
g — > \~ 1 go)- This means that f must become proportional to _1 : f(0 <1)~ -1 f. We see then that f(0) has the 
same asymptotic behaviors as 1//k(0), where /1 is the interpolating function appearing in present MOND theories. 

If q does depend on 0, the EFE does not conform to the standard results. For example, in the DML we could have 
< g(0) ^ 1, in which case scale invariance dictates f(0 < 1) w 0-«(°)f. Then a ~ go{g N ao/ g 2 a ) q{0) = g (v I '°) q(0) ■ We 
see that, as stated above, the condition for external-field dominance, a <C go, when < 1, is indeed always r\ <C 9. 

For example, if g(0) = 1/2, this gives the standard scaling of the MOND acceleration in isolated systems a ~ 
(#»«()) i-e., there is no EFE, except for effects in f (0). So the basic tenets lead to the standard EFE results (indeed 
to an EFE) only if some additional analytic properties are assumed. 

The toy DML theory described by eq.(3), which satisfies scale invariance (but which does not combine with an 
appropriate Newtonian limit), does not lead to an EFE. 

The above analytic assumptions do hold in all the MOND formulations considered to date: e.g., in the original, 
pristine formulation in [1], in the formulation of [3], and in QUMOND [4]. For example, in QUMOND, we can write 
schematically (ignoring the vectorial nature of the quantities involved) 

a/g ~ v [6ii(0) + 9 V M9) + 77] - 1, (16) 



It is difficult to make general statements about the intermediate case where the two accelerations are of the same order. 
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where v{y) is the QUMOND interpolating function, and /j,(x) is such that v[x(i(x)](i(x) = 1. We have /i(9 <C 1) ~ 0, 
fx(0 > 1) w 1; so we see explicitly why the condition 77 <C mm(l,#) is tantamount to a dominant external field. And, 
clearly the next to zcroth-ordcr term is a/ go ~ 77(1 + z>)///(0), where — 1/2 < £> < is the logarithmic derivative of v. 

In summary, the fact that an external field |g | ^> a renders the internal dynamics Newtonian, follows from only 
the basic tenets of MOND, provided that only the instantaneous external field enters the internal dynamics (not 
necessarily true in modified-inertia, time-nonlocal theories). This is relevant, for example, to experimental results in 
the laboratory and solar system, and to the dynamics of star clusters near the sun. On the other hand, the specific 
form of the EFE when |g | <C a , even its very existence, is not strictly dictated by the basic tenets alone. Its basic 
features do follow under another plausible assumption, shared by all full-fledged theories considered to date: that the 
expansion power in eq.(14) does not depend on 9. 

There is no EFE in the DM paradigm. 

G. Disc galaxies have both disc and a halo components of "phantom matter" 

The disc component of a spiral galaxy is described ideally as a thin planar mass. In Newtonian dynamics, the 
surface density of the disc is related to the perpendicular component of the measured acceleration just outside the 
disc, by S = (2-kG)^ 1 g™l (assuming symmetry to reflection in the disc plane). MOND predicts a perpendicular 
component of the acceleration g™l ^ g™±- This difference would be interpreted by a Newtonist as being due to a 
thin disc of 'phantom matter' of surface density 

S = (2^G)- 1 (.g:;l-.g;-), (17) 

just as the general- field difference between the MOND accelerations is interpreted as 'halo of PM'. Such a phantom 
disc is predicted only where there is a baryonic disc, because MOND docs not create a discontinuity when this is not 
forced by the mass distribution. More generally, it would appear to a Newtonist that a certain distribution of disc 
DM is needed to explain the perpendicular (z) dynamics within the disc. 

This effect has been discussed in specific formulations of MOND starting from the first papers on MOND (see, 
in particular, Ref. [29]-where observational indications of disc 'DM' are also discussed-and Ref. [2] with further 
references therein; possible peculiarities of the effect in nonlocal, MI theories were alluded to in Ref. [6]). It is difficult 
to make general statements based only on the basic tenets. However, under some circumstances, the effect can be 
treated similarly to the EFE treated in section IV F above. Consider the perpendicular dynamics up to some height 
Z at galactocentric position r. If the disc were a uniform planar structure, we could have treated the z dynamics 
as a one-dimensional problem with planar symmetry. This might be a good approximation at the very center of the 
disc, provided we explore it to Z <C h, where h is the typical scale length of the disc. In Newtonian dynamics, this 
is a reasonable approximation at all r, provided we limit ourselves to Z <C r. MOND, however does not generally 
permit such a treatment away from the center, even for Z<r: the radial component of the dynamics at r must also 
be reckoned with because of the nonlincarity of MOND. As in the case of the EFE, if we can assume that the only 
effects of the r dynamics on the z dynamics enter through the value of g r (r), the value of the radial acceleration at 
r-which can be assumed a constant as regards the z dynamics at r-then our treatment of the EFE applies here with 
I g 1 replaced by g r (r). 

A disc component of DM, as predicted by MOND, is not a natural occurrence in the cold-dark-matter paradigm. 

H. MOND endows self-gravitating systems with added, but limited stability 

By ND, accelerations in a self gravitating system scale as a ~ pRG (where R is size and p density). In the DML, 
scale invariance (and dimensional considerations) dictates that a ~ (pRflo) 1 ' 2 . So, while in ND Sa/a oc Sp/p, in the 
DML we have 5a/a oc (l/2)5p/p: the response to perturbation is weaker, by a factor ~ 2, leading to higher, but not 
much higher, stability. 

In other words, MOND predicts added stability in DML systems compared with systems of similar structure that 
are in the Newtonian regime. But, it also predicts that the degree of stability does not depend on how deep in the 
MOND regime the systems is. This follows from the invariance of the DML to the two-parameter family of scalings 
discussed in section III. Such scalings relate systems with all possible mass discrepancies, namely with all possible 
values of MG/R 2 ao (as long as they are in the DML), which must then all have the same stability properties. For 
example, all exponential discs, with the same ratio of scale height to scale length, the same ratio of vertical velocity 
dispersion to rotational speeds, etc., but having all possible values of M/R 2 <C a /G, arc predicted by the basic tenets 
of MOND to have the same stability properties. What differs among them is the time scale (for instabilities to grow, 
for example). 
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Detailed discussions of disc stability in specific MOND theories can be found in Refs. [30-34]. Observational 
implications are summarized in Rcf. [2] . 

The predicted MOND dependence of disc stability properties on the mass discrepancy (as evinced by the RC) is not 
reproduced in the DM paradigm: There, a progressively large RC mass discrepancy bespeaks a progressively more 
massive halo compared with the disc. This, in turn, predicts progressively higher disc stability, unless one allows a 
large fraction of the putative DM to reside in a disc component, which is not natural in the cold-dark-matter paradigm. 
In fact, the observations of DML galaxies, with large mass discrepancies, and yet showing signs of instability (such 
as spiral structure and bars) is a challenge for this paradigm. 



I. The excess MOND accelerations (those putatively attributed to DM) cannot much exceed do 

In regions of a system in which g N < a , the excess MOND accelerations 31 cannot much exceed ao- If they 
would, the total acceleration, g, would be much above ao, the dynamics would have to be Newtonian, which would 
contradictorily imply g « g N < ao. This follows essentially from the basic tenets. 32 

However, in systems with g N 3> ao, the maximum value that the excess MOND acceleration can attain depends on 
the fastness with which MOND tends to Newtonian dynamics in the limit ao — > 0. We can demonstrate this with the 
pristine, algebraic MOND relation between g and g N , which can be written as g = v(g N /ao)gN, with v(y Cl)w y -1 ^ 2 , 
v(y > 1) ~ 1. For this case (g — g N )/ao = s {y) = W(y) — l]y, where y = g N /a . So, we are interested in the maximum 
value that s(y) can attain. We see that if v — 1 vanishes slower than a/y (for some constant a), s has no upper 
bound. But, if v— 1 vanishes as a/y or faster, s does have an upper bound. This bound has to be of order one by our 
assumption of weak compatibility. As mentioned in section II, constraints from the solar system force on us 'strong 
compatibility' of MOND with standard dynamics. This dictates that g — g N |/ao vanishes at high accelerations. This 
quantity must thus be bound by a number of order one, which depends on the theory. This prediction was first 
discussed in Ref. [35] in the framework of the nonlinear-Poisson formulation of MOND, and tested and confirmed for 
a sample of disc galaxies in Ref. [36]. 

Let <?^ ax is the maximal Newtonian acceleration in a system. One expects that if g^ ox > ao then the maximal value 
of the excess MOND acceleration would be of order a . But, in a DML system, where g™ ax <C a , the discussion in 
section III implies that the maximum excess acceleration is expected to be g"* ox ~ (aog^"*) 1 / 2 . 

Expressed as a property of a putative DM, this MOND prediction says that the acceleration produced by DM alone 
can never much exceed ao- There is no known reason for such a constraint to hold in the DM paradigm. 



J. Em is a maximum and an accumulation value of the central surface densities of 'dark halos' 

Suppose we have a system of mass density p(r) whose Newtonian acceleration field is g N (r). One then uses some 
techniques to map the dynamics of the system, namely, to measure some components of the accelerations of test 
particles at selected locations in the system, using, e.g., rotation curves, lensing, hot gas hydrostatics, etc. These 
data are then fitted by an acceleration field, forced to be derived from a potential, by making further assumptions 
(for example that the excess acceleration is due to a spherical DM halo of a certain density profile). The MOND 
prediction of this field is g(r) (mind again, MOND does not always predict the existence of an acceleration field; so 
here g is the fit to the MOND prediction). If one attributes the excess g^ = g — g N to DM, then the density of this 
'DM', or what we call 'phantom matter' in the context of MOND, is p p (r) = — (47rG) _1 V • gh- 

Let g max be the maximal characteristic value of \gh\ in the system, and ro the characteristic size of the 'halo': the ra- 
dius at which g/J drops appreciably from g max . Then the characteristic phantom density is (p p ) ~ — (47rG) _1 (V-g;j) ~ 
(4irG)~ 1 g max /ro. So, the characteristic mean surface density of the 'halo' is ~ 2(p p )ro ~ g max /2nG. Seen some- 
what differently, if we assume a quasi-spherical geometry, then p p w — (4:'KG)~ 1 r~ 2 d(r 2 gJ l )/dr = — (AttG)" 1 (dg^/dr + 
2gJ l /r), where g r h is the radial component of g^, (|g[J = |gh|). The two terms are opposite in sign; so if we approximate 



The Newtonian acceleration is always well defined; the 'MOND acceleration' not always so: In MG MOND theories the notion of an 
acceleration field is well defined in the NR case as — V</>, where 4> is the MOND potential dictated by the theory. In MI theories, where 
the actual acceleration of a test particle may depend on details of its orbit, there is no 'acceleration field'. However, as things stand 
today, at least, one always imposes the existence of such a field when interpreting the data. See more details in section IV J. 
Strictly speaking, the second tenet says that dynamics is Newtonian everywhere if g is much larger than ao everywhere. Here I assume 
that this is so also in any sub region in the system. 
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9h/ r ~ —dg^/dr, and integrate we get the same result for (S p ). This relation 

(S p ) ~ ^ maa: S M = ^Sm, (18) 

is a result of the definitions alone. 

Now, according to section IV I, for all systems with g N > a , which is an accumulation value for astrophysical 
systems, we have g ma!C ~ ao, in which case £ ~ 1. 

For DML systems we saw that g ma!C <~ (any™" 21 ) 1 / 2 ; so, 

(S p )=e( 5 r7«o) 1/2 S M , (19) 

with £ ~ 1. The coefficients £, £ ~ 1 depend somewhat on the particular theory (including details, such as the form 
of the relevant interpolating function), and on the exact mass distribution p. However, this leeway in the values of £ 
and £ is small compared with the relative ranges over which system parameters (masses, sizes, etc.) vary. 

This prediction, which concerns a pure 'halo' property, was discussed in detail in Ref. [37], in the context of the 
MG formulation based on the nonlinear Poisson equation, including explicit determination of £ and £. 33 This analysis 
was instigated by the finding of Ref. [39], of a 'universal' halo surface density (which turns out to be very nearly £ M ). 
The ubiquity of this value was already evident in the analysis of Ref. [36] that considered a sample of mostly non 
DML disc galaxies. Additional relevant data and analysis can be found in Ref. [40]. 

A related finding pertaining to the mean baryonic surface density within one 'halo' scale length-which thus constitute 
a relation between the visible and 'dark' matter-was discussed in Ref. [41]. 

If interpreted in the framework of DM, this MOND prediction pertains to pure halo properties. Unlike prediction 
(i), which concerns the asymptotics of the 'halo', the present prediction concerns its bulk properties. Such a critical 
halo surface density is not reproduced in the DM paradigm. 



K. Rotation Curves of Disc Galaxies 



The prediction of RCs of disc galaxies does not count among the 'MOND laws'. Unlike these laws, which each 
pertains to a very limited, one- or two-parameter aspect of galactic dynamics, the RC measures the full midplane 
dynamics (at least for circular orbits). The RC prediction is rather more elaborate, and, in fact, subsumes several of 
the MOND laws: (i), (ii), and (iii). Still, in the present context it is interesting to note that even for RC prediction, 
the room left by the basic tenets for differences between theories is only in the details, not in the essentials. 

This can be seen as follows: two radii characterize the dynamics of a galaxy in MOND: R, the radius containing 
most of the mass, and r M . We saw that the basic tenets fully dictate the RC-in shape and magnitude-far beyond the 
larger of the two (the asymptotic regime) . Now, if R < r M , most of the matter is in the Newtonian regime. The basic 
tenets then dictate a Newtonian-hence theory independent-RC within R, or even far beyond it, if R <C r M (except, 
perhaps, for a small region near the very center). 34 So theories can differ only in the interpolating region just below 
and just above r M . 

When r M < R (strictly speaking r M <C ii), the galaxy is fully in the MOND regime. In this case, characterize 
the disc galaxy, most generally, by a mass M, R, and dimensionless structure parameters r\i (shape parameters, mass 
ratios of different components, etc.). Then, dimensional considerations and scale invariance dictate that the MOND 
rotational speed in the midplane, V(r), can be written, for any MOND theory, in the form 



V(r) = [a rVZ(r)} 1 ' 4 q (r/R,r h ), (20) 

where V N (r) is the Newtonian rotational speed. 35 The dimensionless MG/R 2 ao — (r M /R) 2 cannot appear as an 
argument of the dimensionless and scale-invariant q because it is not scale invariant. Furthermore, q w 1 for r 3> R, 
as dictated by the basic tenets, since the prefactor tends to = (MQq) 1 / 4 there. This is as far as we can take the 
basic tenets on the issue of RCs for DML galaxies. 



From the discussion in Appendix B.2 in Ref. [38], we see that the same results hold in QUMOND (which was used in Ref. [38] without 
recognizing it as QUMOND) 

Depending on the density distribution near the center, the accelerations may drop below ao there. 

Note that while V«, which is defined as Vn = (rgN) 1 ^ 2 , has dimensions of velocity it does not transform like a velocity under scaling, 
since g N scales like A -2 . It is the quantity rV$ that scales appropriately (i.e., is scale invariant). 
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In MI theories, for circular orbits in the midplane of a disc, g = r^ 1 V 2 is necessarily a function of the local value 
g N [5]. Then, for DML discs we have g = (doffw) 1 ^ 2 '? 2 ; this implies that q = 1, irrespective of the mass distribution, 
since the only function of g N with the required properties is a constant (e.g., g N /ao is not scale invariant). 

Some numerical comparisons of DML rotation curves between MI and existing MG theories for several disc models 
were shown in Rcf. [42]. More generally, Ref. [27] showed that the RMS velocities predicted by these two classes in 
the DML differ by no more than five percent. The difference is, however, systematic [V(MI) > V(MG)] and could 
be used to discriminate observationally between the theory classes. 

V. DISCUSSION 

I have listed and discussed a number of predicted 'MOND laws' that should underlie galaxy dynamics. The list 
comprises exact relations, correlations, and constraints that the attributes of galaxies and galactic systems should 
satisfy, according to MOND. These laws are independent in the sense that they do not follow from each other in the 
DM paradigm. Some of them involve ao in various roles, in laws pertaining to baryons alone, to 'DM' alone, or to 
relations between the two components. In all deep-MOND results, ao appears in the combined constant fio = ciqG. 
In other laws, ao appears as a critical surface density £ M = a /27rG, but still in different roles, e.g., as a property of 
baryons in law (iv), and as a property of 'DM' in law (x). 

I showed that to a large extent these laws follow from only the basic tenets of MOND. But there are exceptions, 
where some of the laws, or partial aspects of them, require additional assumptions. Such assumptions are, by and 
large, satisfied by existing MOND theories, but it is not necessary that any MOND theory should satisfy them. 
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